Quantum fluctuations in the effective pseudospin-1/2 model 
for magnetic pyrochlore oxides 
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The effective quantum pseudospin-1/2 model for interacting rare-earth magnetic moments, which 
are locally described with atomic doublets, is studied theoretically for magnetic pyrochlore oxides. 
It is derived microscopically for localized Pr 3+ 4/ moments in VT2TM2O7 (TM =Zr, Sn, Hf, and 
Ir) by starting from the atomic non-Kramers magnetic doublets and performing the strong-coupling 
perturbation expansion of the virtual electron transfer between the Pr 4/ and O 2p electrons. The 
most generic form of the nearest-neighbor anisotropic superexchange pseudospin-1/2 Hamiltonian 
is also constructed from the symmetry properties, which is applicable to Kramers ions Nd 3+ , Sm 3+ , 
and Yb 3+ potentially showing large quantum effects. The effective model is then studied by means 
of a classical mean- field theory and the exact diagonalization on a single tetrahedron and on a 16- 
site cluster. These calculations reveal appreciable quantum fluctuations leading to quantum phase 
transitions to a quadrupolar state as a melting of spin ice for the Pr 3+ case. The model also shows 
a formation of cooperative quadrupole moment and pseudospin chirality on tetrahedrons. A sign of 
a singlet quantum spin ice is also found in a finite region in the space of coupling constants. The 
relevance to the experiments is discussed. 

PACS numbers: 75.10.Jm, 75.10.Dg, 75.30.Kz, 75.50.Ee 
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I. INTRODUCTION 



Quantum fluctuations and geometrical frustration are 
a couple of key ingredients in realizing nontrivial spin- 
disordered ground states without a magnetic dipole long- 
range order (LRO) in three spatial dimensions 1 4 . The 
pyrochlore lattice structure is a typical example where 
the geometrical frustration plays a crucial role in pre- 
venting the LRO 1 ' 5 ' 6 . Of our particular interest in this 
paper is a so-called dipolar spin ice 7 10 , such as Dy 2 Ti 2 07 
and Ho2Ti207, and related systems. The dipolar spin 
ice provides a classical magnetic analogue of a cubic 
water ice 11 and is characterized by the emergent U(l) 
gauge field mediating the Coulomb interaction between 
monopole charges 12 ' 13 as well as the dipolar spin cor- 
relation showing a pinch-point singularity 14 ' 15 . Intro- 
ducing quantum effects to the classical spin ice may 
produce further nontrivial states of matter. Evidence 
of quantum effects has recently been observed with in- 
elastic neutron-scattering experiments on the spin-ice 
related compounds, Tb 2 Ti 2 7 16 18 , Tb 2 Sn 2 7 18 , and 
Pr 2 Sn 2 07 19 . Exploiting a weak coupling of the rare- 
earth magnetic moments to conduction electrons, a chiral 
spin state 3 has been detected through the anomalous Hall 
effect 20 at zero magnetic field without magnetic dipole 
LRO in another related compound Pr 2 Ir 2 07 21 . Vital 
roles of the planar components have also been experimen- 
tally observed in Yb 2 Ti 2 07 and Er 2 Ti 2 07 22 . Obviously, 
quantum fluctuations enrich the otherwise classical prop- 
erties of the spin ice. They may drive it to other states 
of matter, including quadrupolar states and chiral spin 
states 23 . The aims of this paper are to provide a compre- 
hensive derivation of a realistic effective quantum model 
for these spin-ice related materials and to understand its 
basic properties including nontrivial quantum effects. 



A. Classical dipolar spin ice 

Let us briefly review the classical (spin) ice. The low- 
energy properties of water and spin ices are described by 
Ising degrees of freedom that represent whether proton 
displacements (electric dipoles) and magnetic dipoles, re- 
spectively, point inwards ("in") to or outwards ("out") 
from the center of the tetrahedron. The interaction 
among the Ising variables favors nearest-neighbor pairs 
of "in" and "out" and thus suffers from geometrical frus- 
tration. This produces a so-called ice rule 11 ' 24 stabilizing 
"2-in, 2-out" configurations on each tetrahedron. Macro- 
scopic degeneracy of this ice-rule manifold produces Paul- 
ing's residual entropy ^Rln | n . 

In the dipolar spin ice 7 9 , a rare-earth magnetic mo- 
ment rh r = gjfisJr located at a vertex r of tetrahedrons 
plays the role of the Ising variable because of the large 
crystalline electric field (CEF), which is often approxi- 
mately modeled by 
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with the Lande factor gj and £>i s i ng > 0. Here, J is the 
quantum number for the total angular momentum J r , 
and n r defines a unit vector at a pyrochlore-lattice site r 
that points outwards from the center of the tetrahedron 
belonging to one fee sublattice of the diamond lattice and 
inwards to that belonging to the other sublattice. The 
amplitude of the rare-earth magnetic moment is so large 
that the interaction between the magnetic moments is 
dominated by the magnetic dipolar interaction 25 ' 26 , 
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with Ar — r — r' and the summation J2^ r over all the 
pairs of atomic sites. This yields a ferromagnetic cou- 
pling D n .n. = |^ (q/ ^2)3 ~ 2A K between the nearest- 
neighbor magnetic moments for Hk^T^Or and Dy2Ti207 
with the moment amplitude m ~ 10/i# and the lattice 
constant a ~ 10.1 A 9 , providing a main driving force 
of the ice rule. It prevails over the nearest-neighbor su- 
perexchange interaction which is usually assumed to take 
the isotropic Heisenberg form 

n.n. 

#H = — 3J n .n. ^ J r -J r '/J 2 . (3) 

{r,r') 

In the limit of Asing oo, 7f DS i = Rising + H D + 
Hh is reduced to an Ising model 26 , which can explain 
many magnetic properties experimentally observed at 
temperatures well below the crystal-field excitation en- 

er g y 9,10,13,27^ 

Because of the ferromagnetic effective nearest-neighbor 
coupling J e ff = -Dn.n. + Jn.n. > 0, creating "3-in, 1-out" 
and "1-in, 3-out" configurations out of the macroscop- 
ically degenerate "2-in, 2-out" spin-ice manifold costs 
an energy, and can be regarded as defects of mag- 
netic monopoles and anti-monopoles with a unit mag- 
netic charge 6 . Then, the spin ice is described as the 
Coulomb phase of magnetic monopoles where the emer- 
gent U(l) gauge fields mediate the Coulomb interaction 
between monopole charges 12 ' 13 . The density of magnetic 
monopoles is significantly suppressed to lower the total 
free energy at a temperature T < J e ff ~ a few Kelvin. Si- 
multaneously, the reduction of the monopole density sup- 
presses spin-flip processes, for instance, due to a quantum 
tunneling 28 , that change the configuration of monopoles. 
Hence, the relaxation time to reach the thermal equilib- 
rium shows a rapid increase. These phenomena associ- 
ated with a thermal quench of spin ice have been ex- 
perimentally observed 29 ' 30 and successfully mimicked by 
classical Monte-Carlo simulations on the Coulomb gas 
model of magnetic monopoles 27 ' 31 . This indicates that 
the quantum effects are almost negligible in the dipolar 
spin ice. It has been shown that the emergent gapless 
U(l) gauge excitations together with a power-law decay 
of spin correlations can survive against a weak antiferroic 
exchange interaction that exchanges the nearest-neighbor 
pseudospin-1/2 variables ("in" and "out") 12 . This 17(1) 
spin liquid 12 can be viewed as a quantum version of the 
spin ice, though the macroscopic degeneracy of the ice- 
rule manifold should eventually be lifted in the ideal case 
under the equilibrium. 

B. Quantum effects 

At a first glance, one might suspect that quantum fluc- 
tuations should be significantly suppressed by a large to- 
tal angular momentum J of the localized rare-earth mag- 
netic moment and its strong single-spin Ising anisotropy 



Dosing > 0? since they favor a large amplitude of the 
quantum number for J* = J r • n r , either Mj = J 
or —J. Namely, in the effective Hamiltonian, Hbsi — 
Rising + + -#h> a process for successive flips of the 
total angular momentum from Mj = J to — J at one site 
and from —J to J at the adjacent site is considerably 
suppressed at a temperature T <C Di S [ ng . The coupling 
constant for this pseudospin-flip interaction is of order 
of |7 n .n. |(|7 n . n . |/^Di s ing) 2J , and becomes negligibly small 
compared to the Ising coupling J e ff. 

In reality, however, because of the D 3 d crystalline elec- 
tric field (CEF) acting on rare-earth ions [Fig. 1], the 
conservation of J 2 , which is implicitly assumed in the 
above consideration, no longer holds in the atomic level. 
Eigenstates of the atomic Hamiltonian including the LS 
coupling and the CEF take the form of a superposition 
of eigenstates of J z whose eigenvalues are different by in- 
teger multiples of three. Obviously, this is advantageous 
for the quantum spin exchange to efficiently work. 

Attempts to include quantum effects have recently 
been in progress. It has been argued that the presence 
of a low-energy crystal-field excited doublet above the 
ground- state doublet in Tb2Ti207 10 ' 16 enhances quan- 
tum fluctuations and possibly drives the classical spin 
ice into a quantum spin ice composed of a quantum su- 
perposition of "2-in, 2-out" configurations 32 . 

We have proposed theoretically an alternative scenario, 
namely, a quantum melting of spin ice 23 . A quantum en- 
tanglement among the degenerate states lifts the macro- 
scopic degeneracy, suppresses the spin-ice freezing, and 
thus leads to another distinct ground state. Actually, the 
quantum-mechanical spin-exchange Hamiltonian mixes 
"2-in, 2-out" configurations with "3-in, 1-out" and "1- 
in, 3-out", leading to a failure of the strict ice rule and 
a finite density of monopoles at the quantum-mechanical 
ground state. Namely, the quantum-mechanically prolif- 
erated monopoles can modify the dipolar spin-ice ground 
state, while a spatial profile of short-range spin cor- 
relations still resembles that of the dipolar spin ice 23 . 
They may appear in bound pairs or in condensates. We 
have reported that there appears a significantly large 
anisotropic quantum-mechanical superexchange interac- 
tion between Pr magnetic moments in Pr2TM207 23 (TM 
= Zr, Sn, Hf, and Ir) 33 . This anisotropic superexchange 
interaction drives quantum phase transitions among the 
spin ice, quadrupolar states having nontrivial chirality 
correlations, and the quantum spin ice, as we will see 
later. 

Actually, among the rare-earth ions available for mag- 
netic pyrochlore oxides 10 ' 33 , the Pr 3+ ion could optimally 
exhibit the quantum effects because of the following two 
facts, (i) A relatively small magnitude of the Pr 3+ lo- 
calized magnetic moment, whose atomic value is given 
by 3.2/iB, suppresses the magnetic dipolar interaction, 
which is proportional to the square of the moment size. 
Then, for Pr 2 TM 2 7 , one obtains D> n .n. ~ 0.1 K, which 
is an order of magnitude smaller than 2.4 K for Ho2Ti207 
and Dy2Ti207. Similarly, quantum effects might appear 
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prominently also for Nd 3+ , Sm 3+ , and Yb 3+ ions be- 
cause of their small moment amplitudes, 3.3/i^, 0.7/zb, 
and 4//#, respectively, for isolated cases, (ii) With fewer 
4/ electrons, the 4/-electron wavefunction becomes less 
localized at atomic sites. This enhances the overlap with 
the O 2p orbitals at the 01 site [Fig. 1 (a)], and thus 
the superexchange interaction which is also further in- 
creased by a near resonance of Pr 4/ and O 2p levels. 
Moreover, this superexchange interaction appreciably de- 
viates from the isotropic Heisenberg form because of the 
highly anisotropic orbital shape of the /-electron wave- 
function and the strong LS coupling. Since the direct 
Coulomb exchange interaction is even negligibly small 25 , 
this superexchange interaction due to virtual f-p electron 
transfers is expected to be the leading interaction. 

Recent experiments on P^S^Or 34 , P^Z^Ot 35 , and 
Pr2lr2 07 36 have shown that the Pr 3+ ion provides the 
(111) Ising moment described by a non-Kramers mag- 
netic doublet. They show similarities to the dipolar 
spin ice. (i) No magnetic dipole LRO is observed down 
to a partial spin- freezing temperature Tf ~ 0.1-0.3 
K 19 ' 21 ' 34 37 . (ii) Pr2lr207 shows a metamagnetic tran- 
sition at low temperatures only when the magnetic field 
is applied in the [111] direction 21 , indicating the ice- 
rule formation due to the effective ferromagnetic cou- 
pling 2 J e ff ~ 2J n . n . ~ 1.4 K 21 . On the other hand, sub- 
stantially different experimental observations from the 
dipolar spin ice have also been made. The Curie- Weiss 
temperature Tew is antiferromagnetic for the zirconate 35 
and iridate 36 , unlike the spin ice. The stannate shows a 
significant level of low-energy short-range spin dynamics 
in the energy range up to a few Kelvin 19 , which is absent 
in the classical spin ice. Furthermore, the iridate shows 
the Hall effect at zero magnetic field without magnetic 
dipole LRO 21 , suggesting an onset of a chiral spin-liquid 
phase 3 at Tu ^ 1.5 K. 

The discovery of this chiral spin state endowed with 
a broken time-reversal symmetry on a macroscopic scale 
in Pr2lr207 without apparent magnetic LRO 21 has in- 
creased the variety of spin liquids. One might speculate 
that this is caused mainly by a Kondo coupling to Ir 
conduction electrons and thus the RKKY interaction 38 . 
However, the low-temperature thermodynamic proper- 
ties are common in this series of materials, Pr2TM207, 
except that a small partial reduction (~ 10%) of Pr mag- 
netic moments probably due to conduction electrons af- 
fects the resistivity and the magnetic susceptibility in 
Pr2lr2 07 36 . Furthermore, the onset temperature Th ~ 
1.5 K for the emergent anomalous Hall effect is compa- 
rable to the ferromagnetic coupling 2J e ff ~ 1.4 K 21 . 
Therefore, it is natural to expect that a seed of the chiral 
spin state below Tjj exists in the Pr moments interacting 
through the superexchange interaction and possibly the 
state is stabilized by the conduction electrons. Another 
intriguing observation here is that without appreciable 
quantum effects, the chiral manifold of classical ice-rule 
spin configurations 21 that has been invented to account 
for the emergent anomalous Hall effect will result from 



a magnetic dipole LRO or freezing, which is actually ab- 
sent down to Tf. This points to a significant level of 
nontrivial quantum effects. 




FIG. 1: (Color online) (a) Pr 3+ ions form tetrahedrons 
(dashed lines) centered at O 2- ions (01), and are sur- 
rounded by O 2- ions (02) in the D 3 d symmetry as well as 
by transition-metal ions (TM). Each Pr magnetic moment 
(bold arrow) points to either of the two neighboring 01 sites. 
(x r ,y r , z r ) denotes the local coordinate frame, (b) The local 
coordinate frame from the top. The upward and the down- 
ward triangles of the O 2- ions (02) are located above and 
below the hexagon of the TM ions. 

In this paper, we develop a realistic effective theory 
for frustrated magnets Pr 2 TM 2 07 on the pyrochlore lat- 
tice and provide generic implications on quantum effects 
in spin-ice related materials, giving a comprehensive ex- 
planation of our recent Letter 23 . In Sec. II, the most 
generic nearest-neighbor pseudospin-1/2 Hamiltonian for 
interacting magnetic moments on the pyrochlore lattice 
is derived on a basis of atomic magnetic doublets for both 
non-Kramers and Kramers ions. In particular, it is mi- 
croscopically derived from strong-coupling perturbation 
theory in the Pr 3+ case. We analyze the model for the 
non-Kramers case by means of a classical mean-field the- 
ory in Sec. Ill, which reveals spin-ice, antiferroquadrupo- 
lar, and noncoplanar ferroquadrupolar phases at low tem- 
peratures. Then, we perform exact-diagonalization cal- 
culations for the quantum pseudospin-1/2 case on a single 
tetrahedron in Sec. IV and on the 16-site cube in Sec. V. 
We have found within the 16-site cluster calculations a 
cooperative ferroquadrupolar phase, which is accompa- 
nied by crystal symmetry lowering from cubic to tetrag- 
onal and can then be categorized into a magnetic analog 
of a smectic or crystalline phase 39 . This provides a sce- 
nario of the quantum melting of spin ice and can explain 
the experimentally observed magnetic properties, includ- 
ing powder neutron-scattering experiments on Pr 2 Sn 2 07 
and the magnetization curve on Pr2lr207. We also reveal 
a possible source of the time-reversal symmetry breaking 
observed in Pr2lr207. It takes the form of the solid angle 
subtended by four pseudospins on a tetrahedron, each of 
which is composed of the Ising dipole magnetic moment 
and the planar atomic quadrupole moment, and shows 
a nontrivial correlation because of a geometrical frus- 
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tration associated with the fee sublattice structure. A 
possible sign of a singlet quantum spin-ice state has also 
been obtained within the 16-site numerical calculations 
in another finite region of the phase diagram. Sec. VI is 
devoted to discussions and the summary. 

II. DERIVATION OF THE EFFECTIVE MODEL 

In this section, we will give a microscopic derivation 
of the effective pseudospin-1/2 Hamiltonian in a com- 
prehensive manner. Though we focus on 4/ localized 
moments of Pr 3+ ions, the form of our nearest-neighbor 
anisotropic pseudospin-1/2 Hamiltonian is most generic 
for atomic non-Kramers magnetic doublets. We will also 
present the generic form of the nearest-neighbor Hamil- 
tonian for Kramers doublets of Nd 3+ , Sm 3+ , and Yb 3+ . 

A. Atomic Hamiltonian for Pr 3+ 

1. Coulomb repulsion 

The largest energy scale of the problem should be the 
local Coulomb repulsion among Pr 4/ electrons. A pho- 
toemission spectroscopy on Pr 2 TM 2 07, which is desir- 
able for its reliable estimate, is not available yet. A typi- 
cal value obtained from Slater integrals for Pr 3+ ions is of 
order of 3-5 eV 40 . A cost of the Coulomb energy becomes 
0, U, and 3U for the occupation of one, two, and three 
/ electrons, respectively [Fig. 2). For Pr 3+ ions, the O 
2p electron level A at the 01 site should be higher than 
the f 1 level and lower than the f 3 level [Fig. 2]. Then, it 
would be a reasonably good approximation to start from 
localized f 2 states for Pr 3+ configurations and then to 
treat the other effects as perturbations. 

f 
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FIG. 2: (Color online) Local level scheme for / and p elec- 
trons, and the local quantization axes z r and z r /. 



2. LS coupling for f 2 configurations 

We introduce operators J, £, and S for the total, or- 
bital, and spin angular momenta of f 2 electron states 
of Pr 3+ . Within this f 2 manifold, the predominant LS 
coupling Xls > in Hls — ^lsL ■ S gives the ground- 
state manifold 3 #4 with the quantum numbers J = 4, 
L = 5, and S = 1 for the total, orbital, and spin angular 
momenta, respectively. 



3. Crystalline electric field 

The ninefold degeneracy of the ground-state manifold 
3 i^4 is partially lifted by the local crystalline electric field 
(CEF), which has the Dsd symmetry about the (111) di- 
rection toward the 01 site. We define the local quanti- 
zation axis z r as this (111) direction. Then, the Hamil- 
tonian for the CEF, 

3 

#CEF = ^2 V CEF l1 ^2 ^2 f ] r,m u afr,m' v a, (4) 
mi,m f l =— 3 r a=± 

contains not only orthogonal components with mi = m\ 
but also off-diagonal components with raj — raj = ±3 and 
±6, all of which become real if we take x and y axes as x r 
and y r depicted in Figs. 1 (a) and (b). Here, f r mi a and 
fr mi a denote the annihilation and creation operators of 
an / electron with the z components mi and m s = a/2 
of the orbital and spin angular momenta, respectively, in 
the local coordinate frame at the Pr site r. The formal 
expressions for V^^™ 1 within the point-charge analysis 
are given in Appendix A. In the rest of Sec. II A, we drop 
the subscript for the site r for brevity. 

We perform the first-order degenerate perturbation 
theory, which replaces Eq. (4) with P( 3 H a )H C efP{ 3 H 4 ), 
where P^H^) is the projection operator onto the 
3 i^4 manifold. First, let us introduce a notation of 
|L, Ml] S, Ms) for the f 2 eigenstate corresponding to the 
orbital and spin quantum numbers (L,Ml) and (S,Ms) 
in the local coordinate frame. It is straightforward to ex- 
press the eigenstates {\Mj)}mj=-j,— ,j of J z in terms of 
1 5, Ml; 1, Ms) and then in terms of /-electron operators, 

\Mj) = Ml Ms 

M L Ms 

M L Ms 

as explicitly written in Appendix B. 

Finally, we obtain the following representation of 
Eq. (4) in terms of Eq. (5) 
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3 

(Mj\H C Ef\Mj) = ^ Vcef" 1 ^ [CMj.mi.ml'.ff.ff'C'Ml.mj.mj'.ff.ff' - C'Mj.mj'.mj.ff.ff'C'. 



mi ,m( ,m^ =— 3 



m" ,mi ,a,a ,{ ^ M'j^m^ ,m" ,0-' ,cr 



(6) 



The CEF favors Mj = ±4 configurations that are lin- 
early coupled to Mj = ±1 and ^2 because of the Dsd 
CEF. This leads to the atomic non-Kramers magnetic 
ground-state doublet, 

\a) D = a\4a)+pa\a)-j\-2a), (7) 

with small real coefficients (3 and 7 as well as a = 
y/l - f3 2 -7 2 . For Pr 2 Ir 2 7 , the first CEF excited state 
is a singlet located at 168 K and the second is a doublet 
at 648 K 41 . They are located at 210 K and 430 K for 
Pr 2 Sn 2 07 19 . These energy scales are two orders of mag- 
nitude larger than that of our interest, 2J n , n ~ 1.4 K. 
Hence it is safe to neglect these CEF excitations for our 
purpose. Then, it is convenient to introduce the Pauli 
matrix vector fr r for the pseudospin-1/2 representing the 
local doublet at each site r, so that Eq. (7) is the eigen- 
state of cr* = & r • z r with the eigenvalue a. 

Note that in the case of Tb 3+ , the first CEF excited 
state is a doublet at a rather low energy ~ 18.7 K 42 , 
and the effects of the first excited doublet cannot be ig- 
nored 32 when a similar analysis is performed. Neverthe- 
less, since this CEF excitation in Tb 3+ is an order of 
magnitude larger than J e ff, it could be integrated out 32 . 
Then, the model reduces to a similar form of the effec- 
tive pseudospin-1/2 Hamiltonian which has been derived 
in Ref. 23 and will also be discussed below, though the 
explicit form has not been presented as far as we know. 

B. Dipole and quadrupole moments 

Only the z component a* — & r • z r of the pseudospin 
contributes to the magnetic dipole moment represented 
as either "in" or "out" , while the transverse components 
a x — & r • x r and = cr r • y r correspond to the atomic 
quadrupole moment, i.e., the orbital. This can be eas- 
ily shown by directly calculating the Pr 4/ magnetic 
dipole and quadrupole moments in terms of the pseu- 
dospin. We first take the projection of the total angular 
momentum J = ( J x , J y , J z ) to the subspace of the local 
non-Kramers magnetic ground-state doublet described 
by Eq. (7). It yields 

d (o-\J z \ct') d = (4a 2 +/3 2 - 2 7 2 )a^, (8) 
dHJV^ = 0, (9) 

with J ± = (J x ± ijy). With the Lande factor gj = 4/5 
and the Bohr magneton the atomic magnetic dipole 



moment is given by 

rhr = gjl^B^a 2 +f3 2 - 2 7 2 )<j> r . (10) 

Note that cannot linearly couple to neutron spins 
without resorting to higher CEF levels. On the other 
hand, the quadrupole moments are given by 

D (a\{J z J ± }\a , ) D = -3607^,-*'. (H) 

This is a general consequence of the so-called non- 
Kramers magnetic doublet and not restricted to the Pr 3+ 
ion. Namely, when the atomic ground states of the 
non-Kramers ions having an even number of / electrons 
and thus an integer total angular momentum J are de- 
scribed by a magnetic doublet, only a z r contributes to the 
magnetic dipole moment, while cr x,y corresponds to the 
atomic quadrupole moment. This sharply contrasts to 
the following two cases: (i) In the case of Kramers dou- 
blets, all three components of & r may contribute to the 
magnetic dipole moment while their coefficients can be 
anisotropic, (ii) In the case of non-Kramers non-magnetic 
doublets, & r corresponds to a quadrupole moment or 
even a higher-order multipole moment which is a time- 
reversal invariant. 



C. Superexchange interaction 

Now we derive the superexchange Hamiltonian through 
the fourth-order strong-coupling perturbation theory. 
Keeping in mind the local level scheme of Pr 4/ electrons 
and O 2p electrons at 01 sites, which has been explained 
in Sec. II A, we consider nonlocal effects introduced by 
the electron transfer between the Pr 4/ orbital and the 
O 2p orbital. 

1 . Local coordinate frames 

In order to symmetrize the final effective Hamiltonian, 
it is convenient to choose a set of local coordinate frames 
so that it is invariant under 180° rotations of the whole 
system about three axes that include an 01 site and are 
parallel to the global X, F, or Z axes, which belong to 
the space group F d ^ m of the present pyrochlore system. 
We can start from the local coordinate frame previously 
defined in Sec. II A 3 and in Fig. 1 for a certain site and 
generate the other three local frames by applying the 
above three rotations. 
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For instance, we can adopt 

*o = -^(1,1,-2), 

y = ^ (-1,1,0), 

z = ^(M,l), (12a) 
for the Pr sites at R + ao with ao = — f (1, 1, 1), 
xi = -^(1,-1,2), 

Vi = ^ (-1,-1,0), 

*i = ^(1,-1,-1), (12b) 
for the Pr sites at R + a\ with a\ = | (— 1, 1, 1), 
*2 = (-1,1,2), 
y 2 = ^ (1,1,0) , 

*2 = ^(-1,1,-1), (12c) 
for the Pr sites at R + a,2 with a2 = f (1, —1, 1), and 
x 3 = -±=(-1,-1,-2), 

y 3 = ^(i,-i,o), 

zs = -^(-1,-1,1), (12d) 

for the Pr sites at R + as with a3 = |(1, 1, — 1), where 
R represents a fee lattice vector R = J2i=i 2 3 n i^i 
spanned by R\ = (0, a/2, a/2), = (a/2,0, a/2), and 
R3 = (a/2, a/2, 0) with integers (ni, ri2, 723), and a is the 
lattice constant, i.e., the side length of the unit cube. In 
particular, all the local z directions attached to the Pr 
sites belonging to the tetrahedron centered at the 01 site 



R point inwards, and they satisfy the relation 

3 

J2(*i,yuZi) = (0,0,0). (13) 

i=0 

Actually, other sets of local coordinate frames which are 
obtained by threefold and sixfold rotations about [111] 
yield exactly the same expression for the effective Hamil- 
tonian for Kramers and non-Kramers cases, respectively. 

These local coordinate frames are related to the fol- 
lowing rotations of the global coordinate frame, 

R( r \<Pi,#i) = ('a^W**), (14a) 
ipo = 7r/4, $0 = arccos , (14b) 

(f! = 3tt/4, #1 = -7T + arccos (l/V^j , (14c) 

(p 2 = -7r/4, $ 2 = -7r + arccos (l/V^j , (14d) 

(fs = -3tt/4, tf 3 = arccos (l/V^j . (14e) 

Note that the coordinate frame for the spins is always 
attached to that for the orbital space in each case. The 
rotation of j = I + s with the orbital I and the spin s of 
a single electron takes the form 



exp 



-itpi] exp 



(15) 



2. f-p hybridization 

The 4/ electrons occupying the atomic ground-state 
doublet, Eq. (7), or 4/ holes can hop to the O 2p levels 
at the neighboring 01 site. Because of the symmetry, 
the f-p electron transfer along the local z axis is allowed 
only for the pfa bonding (mi = 0) between / 2 (5^ 2 -3r 2 ) 
and p z orbit als and the pfn bondings (mi = ±1) between 
f x {bz 2 -r 2 ) an d p x orbitals and between f y (5 Z 2- r 2) and p y 
orbitals in the local coordinate frame defined in Eqs. (12) 
[Fig. 3 (a)]. Their amplitudes are given by two Slater- 
Koster parameters 43 V p f a and V p fn, respectively. Then, 
the Hamiltonian for the f-p hybridization reads 



Refccr=± m z ,raj=0,±l 



<7,<r'=± 



ai,mi,a \ ^R-\-ai) PR+{1+t)cl 

\ y mi ,m, ;a,(T 



(16) 



with V±i = V^/tt, Vo = V p f a , and an index r = ± for 
the two fee sublattices of the diamond lattice, where 
PR+(i+ T )ai,mi,a represents the annihilation operator of 
a 2p electron at the 01 site R + (1 + r)ai with the or- 



bital and spin quantum numbers m\ and m s = a/2, re- 
spectively, in the global coordinate frame. Here, Rji+ a . 
transforms the representation from the global frame for 
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(a) 




FIG. 3: (Color online) (a) Two /-p transfer integrals; V p f a 
between p z and f^ z 2_ 3r 2^ z orbitals, and V v f^ between p x /p y 
and f x (5 Z 2 _ r 2^/ f y ( 5z 2_ r 2y (b) f-p virtual electron hopping 

processes, n (n) and t in the state f n p £ f n represent the 
number of / electrons at the Pr site r (r r ) and that of p 
electrons at the 01 site. 



3. Strong- coupling perturbation theory 



Now we are ready to perform the strong-coupling per- 
turbation expansion in V p f n and V p f a . Hybridization 
between these Pr 4/ electrons and O 2p electrons at the 
01 site, which is located at the center of the tetrahedron, 
couples f 2 states having the local energy U with f 1 and 
f 3 states having the local energy levels and 3/7, respec- 
tively [Fig. 2]. Here, the LS coupling has been ignored 
in comparison with U for simplicity. Creating a virtual 
p hole decreases the total energy by A, which is the p 
electron level measured from the f 1 level. 

First, the second-order perturbation in V p f a and V p f n 
produces only local terms. They only modify the CEF 
from the result of the point-charge analysis with renor- 
malized parameters for the effective ionic charges and 
radii. Nontrivial effects appear in the fourth order in 
V p f a and Vpfn- Taking into account the virtual processes 
shown in Fig. 3 (b), the fourth-order perturbed Hamilto- 
nian in V p f a and V p f^ is obtained as 



H ff = 



(orr _ v mi v m f i v m2 v m f 2 p r:rniiCri f^ 

{r,r')rni,rn2rn' 1 ,rn' 2 cri,<J2<j l 1 ,<j' 2 



2U 



-^S mum2 5 m ^ ml2 5 aua2 5 ai y 2 + ( 2 jj_ A + Jj) (4^) miimi; „ liff ; ( R r' Rr )m' 1 



m2',cr 1 ,<J2 



(17) 



D. Effective pseudospin-1/2 model 

Next we project the superexchange Hamiltonian 
Eq. (17) onto the subspace of doublets given by Eq. (7). 
For this purpose, we have only to calculate for a site r 
the matrix elements of the operators P r ,mi,afr,m' ,a' with 
m/,mj = 0,±1 and <r, a' — ±1, in terms of \Mj) that is 
explicitly represented with /-electron operators in Ap- 
pendix B, and then in terms of the atomic doublet \<j)d, 
Eq. (7). Then, we finally obtain the effective quantum 
pseudospin- 1/2 Hamiltonian; 

n.n. 

#eff = Jn.n. ^ [KK' + 25 (°r K' + °r °r>) 
(r,r f ) 

+2g(e 2 ^(j+a++/i.c.)], (18) 

with &r = (a* ± i&y)/2, where fr r represents a vector of 
the Pauli matrices for the pseudospin at a site r. The 
phase 44 (j) r y takes — 27r/3, 27r/3, and for the bonds 
shown in blue, red, and green colors in Fig. 4, in the local 
coordinate frames defined in Eq. (12). This phase can- 
not be fully gauged away, because of the noncollinearity 
of the (111) magnetic moment directions and the three- 
fold rotational invariance of (r, cr r ) about the [111] axes. 



Equation (18) gives the most generic nearest-neighbor 
pseudospin-1/2 Hamiltonian for non-Kramers magnetic 
doublets of rare-earth ions such as Pr 3+ and Tb 3+ that is 
allowed by the symmetry of the pyrochlore system. Note 
that the bilinear coupling terms of a% and a^, are prohib- 
ited by the non-Kramers nature of the moment; namely 
a z r changes the sign under the time-reversal operation, 
while a%' y does not. 

The dependence of the Ising coupling constant on C7, 
A, V p f a and V p f n takes the form, 

J n.n. = (2/ /_" A)2 (jj + 2/7- A J 7 ' V Pf*/ V Pf^> 

(19) 

where J contains the dependence on the remaining di- 
mensionless variables 7, V v f^jVpfa)- We show the 
sign of J as functions of 7 and Vpf^/Vpfa for several 
choices of /3/j in Fig. 5. In particular, for —0.37 < 
Vpfjr/Vpfa < —0.02 which includes a realistic case of 
Vpf n /V p fa ~ —0.3, J is found to be positive. Since the 
prefactor in Eq. (19) is positive, the Ising coupling J n . n . is 
also positive, namely, antiferroic for pseudospins, in this 
case. Taking account of the tilting of the two neighboring 
local z axes by = arccos(— 1/3), this indicates a ferro- 
magnetic coupling between the physical (111) magnetic 
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the cases of = 0, ±1/3, ±1/6. Henceforth, we adopt 
/3 = 7.5% and = 1/3, following the point-charge 
analysis of the inelastic powder neutron-scattering data 
on Pr2lr207 [Ref. 41]. Actually, these estimates of j3 
and 7 lead to the local moment amplitude, 

M = gjfis^a 2 ± /3 2 - 2 7 2 ) « 2.9/xb, (20) 

according to Eq. (10), which reasonably agrees with the 
experimental observation on P^I^Or 36 and P^Z^Or 35 . 
Then, we obtain S ~ 0.51 and q ~ 0.89, indicating the 
appreciable quantum nature. In general, however, the 
values of S and q may vary depending on a transition- 
metal ion and crystal parameters. We note that a finite 
q was not taken into account in the literature seriously. 
It is instructive to rewrite the q term as 

q ((<7 r • n r y){v r > ■ n r y) - (a r • n' r r ,)(a r • raj.,r')) -( 21 ) 

where we have introduced a two-dimensional vector com- 
posed of the planar components of the pseudospin, B r = 
(a^,a^), and two orthonormal vectors 



FIG. 4: (Color online) The pyrochlore lattice structure. The 
phase 4> r ,r' appearing in Eq. (18) takes — 27r/3, 27r/3, and 
on the blue, red, green bonds, respectively, in our choice of 
the local coordinate frames [Eqs. 12]. 



P=-y/3 

(3=-y/6 

p = o 

p = y/6 

(3 = y/3 



-0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 

Y 

FIG. 5: (Color online) A phase diagram for the sign of the 
dimensionless Ising coupling J, defined through Eq. (19), as 
functions of 7 and Vpf^/Vpfa for ft = 7/3, 7/6, 0, —7/6, and 
—7/3. In each case, J is positive in the shaded region, and 
negative otherwise. 



fi r , r ' — (COS 4> r ,r', ~ Sm r?r / ) , (22a) 

n' r r , = (sm(j) r ^ r ' ,cos(j) r y). (22b) 

Then, it is clear that the sign of q can be absorbed by 
rotating all the pseudospins cr r about the local z r axes 
by 7r/2. Furthermore, in the particular case 5 — \q\ or 
— the planar (<j r ) part, namely, the sum of the S and 
q terms, of the Hamiltonian Eq. (18) is reduced to the 
antiferroic or ferroic pseudospin 120° Hamiltonian 45 . 



— - - - - (3 = y/3 (b) q 

P = y/6 I 

P=0 \ ' / 




-0.2 



moments and provides a source of the ice rule. 

The D^d CEF produces two quantum-mechanical 
interactions in the case of non-Kramers ions; the 
pseudospin-exchange and pseudospin-nonconserving 
terms. The ratios S and q of their coupling constants 
to the Ising one are insensitive to U/V p f a and A/V p f a 
but strongly depends on f3 and 7. Figures 6 (a) and (b) 
show 5 and q, respectively, as functions of 7 that charac- 
terizes the Dsd CEF for a typical choice of parameters, 
u / v pfcr = 5, A/Vpf a = 4, and V pf7r /V p f a = -0.3, in 



FIG. 6: (Color online) The coupling constants (a) 8 and (b) 
q as functions of 7 for several choices of = 1/3, 1/6, 0, 
-1/6, and -1/3. We have adopted U/V pfa = 5, A/V pftT = 4, 
and Vpfn/Vpfa -0.3. 

For Kramers ions such as Nd 3+ , Er 3+ , and Yb 3+ , there 
appears another coupling constant 46,47 for an additional 
interaction term 

n.n. 

H K = K ^2 [°r (^V * n ri r f> j ± (v r ' ^r.r^j K> : ( 23 ) 
(r,r') 
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whose form has been obtained so that it satisfies the 
threefold rotational symmetry about (111) axes, i.e., z r , 
the mirror symmetry about the planes spanned by z r and 
z r i for all the pairs of nearest- neighbor sites r and r', 
and the twofold rotational symmetry about X, Y, and 
Z axes. This reflects the fact that all the components of 
& r change the sign under the time-reversal, & r — » — & ri 
and hence Hk respects the time-reversal symmetry for 
Kramers ions. Equation (18) and the additional term 
Eq. (23) appearing only for Kramers doublets define the 
most generic form of the nearest-neighbor bilinear inter- 
acting pseudospin-1/2 Hamiltonian that is allowed by the 
symmetry of the magnetic pyrochlore system. Through- 
out this paper, we restrict ourselves to the case of non- 
Kramers ions, for which K vanishes and Hk does not 
appear. 



III. CLASSICAL MEAN-FIELD THEORY 

We start with a classical mean-field analysis of our ef- 
fective Hamiltonian Eq. (18) along the strategy of Ref. 5. 
We look for the instability with decreasing temperature, 
and then consider candidates to the mean-field ground 
state in the space of the two coupling constants S and 
q. We restrict ourselves to the nearest-neighbor model, 
though it is known that longer-range interactions can lift 
the degeneracy at least partially 5 . Since a* and a^r are 
decoupled in the Hamiltonian in the classical level, we 
proceed by requiring that either of (a*) and/or (a% ,y ) 
is finite. It reveals three distinct mean-field instabili- 
ties in the case of J n . n . > 0, which also provide candidate 
mean-field ground states under the constraint |(<x r )| < 1. 
Note that the decoupling approximation of the Ising and 
planar components becomes inaccurate when the SU(2)- 
symmetric point of Heisenberg antiferromagnet. 



A. Ising states (<r£) ^ 
Let us introduce a vector, 



R+clq J ^R+ai > u R+a 2 ' u R+a 3 



>*fn-«s)> ( 24 ) 



where R is a fee lattice vector and {a^}^o, - ,3 have been 
defined in Sec. II CI. In the mean-field approximation, 



magnetic dipolar states characterized by a nonzero (cr*) 
is obtained as the states having the minimum eigenvalue 
of the following mean-field Hamiltonian, 



n z MF = N T J2{dl)h z q (d q ) 



(25) 



fq y +q z fq z +q x fq x +q y 

2J nT [ fq x -q y fq z -q x l^g) 

fq z +q x fq x ~q y fq y ~q z 



where f q = cos(ga/4), and 



q x +q y fq x ~q 



fq 







/ <t£. p -iq-{R+a ) \ 
/ a R+a e \ 



Nn 



R 



J R+ax 



iq-(R+a ) 
e iq-(R+a!) 
p -iq-(R+a 2 ) 



(27) 



with Nt = N/4 where TV is the number of pyrochlore- 
lattice sites. h z has the eigenvalues 5 , 



e\ = -2 J n . n . , 2 J n . n . (1 =f y/1 + 9q) 



(28) 



Where g q = f 2 qj2q y + f2q y f2q z + f2qj2q x . For Jn.n. > 

0, the lowest energy of the mean-field solution, which is 
obtained as — 2J n . n . per tetrahedron for any wavevector 
q all over the Brillouin zone reflecting the macroscopic 
degeneracy, coincides with the exact ground-state energy 
of the nearest-neighbor spin ice model. 



B. Planar states (a=) / 

Introducing another vector, 

Qr = (<?R+ ao i CR+ai i &R-\-cl2 i &R+CL3 ), (29) 

with a r = the mean-field Hamiltonian reads 

^MF = £<<%>^<&>, (30) 

q 

where 



-27 





-VSr x -r z \ £ 

2 )Jq y +q z 







(St + q- 

(Sf + q^=^)f qi+q , (Sf + qt z )f qx . qy 
V (Sr + qr z )f qx+qy (6r + f. 



{S f + q ^^)f qz+qx (5r + qr z )f qx+qy 
(Sfo + qf z )U- qy (5f + £d^ )Uz _ qx 
{5ro+ -^ f r^ )f qy - qz 



0/< 



q y -q z 



2 





r 



We have introduced the Fourier component Q q of Qr in For S > (\q\ + l)/2, the lowest eigenvalue of is given 



an analogy to Eq. (27). 



by 



_PAF 



-2J n . n .(<5 + 2|g|) 



(32) 
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for the planar antiferro-pseudospin (PAF) states at the 
rods q = ^^(l, ±1, ±1) with h being an arbitrary real 
number. It has the 120° planar pseudospin structure 
within a tetrahedron, which is expressed as 



f °> < ,ai > H a ,a 2 , H a ,a 3 ) ^ Q = Qx(h 1, 1) 
I < ,ai > ^ n'a u a 2 ,n' ai ,a 3 ) for q = ^(-1, 1, 1) 



for gr = ^(1,-1,1) 



(33) 

for q > 0. When g < 0, n' in Eq. (33) are replaced by n. 

In the other case of 5 < (\q\ + l)/2, the lowest eigen- 
value of hq is given by 



(34) 



for the planar ferro-pseudospin (PF) state at q = and 
symmetry-related q vectors connected by (^i, ^2,^3) 
and/or ^(1, 1, 1) with integers 774, 712, and 713. This state 
has eigenvectors showing a collinear ferroic alignment of 
the planar components of the pseudospins. 

Because of the saturation of each ordered moment, ex- 
cept at one site per tetrahedron in the case of the PAF 
phase, these states can be stabilized as the ground state 
if the energy is lower than in the dipolar state. 



\q\ 



PF (AFQ) 
Q=0 



PAF Qrod//[111] 

(noncoplanar FQ for q>0) 
(AFQforq>0) 



Ising spin Ice 



-1/3 



A N.N. spin ice 

* N.N. Heisenberg AF 



FIG. 7: (Color online) Classical mean-field phase diagram 
of the model given by Eq. (18). The Q = planar ferro- 
pseudospin (PF) phase physically represents the antiferro- 
quadrupole (AFQ) phase. The planar antiferro-pseudospin 
(PAF) phase is characterized by Bragg rods Q || [111] and 
physically represents a 2D AFQ phase for q > or a 2D 
noncoplanar ferroquadrupole (FQ) phase for q < 0, where 
quadrupole moments are aligned only within the plane per- 
pendicular to Brag rod vectors Q || [111] in the mean-field 
level. In the region around the N.N. Heisenberg antiferro- 
magnet, the present mean-field theory becomes less accurate. 
The point X = (6, q) = (0.51, 0.89) for Pr 3+ is also shown. 



C. Mean-field phase diagram 

Comparing the energies of the dipolar state and the 
two quadrupolar states given above, we obtain the clas- 
sical mean-field phase diagram shown in Fig. 7; (i) 
the macroscopically degenerate dipolar states associated 
with the nearest-neighbor spin ice for —1/3 < S < 1— 2\q\, 
(ii) the planar antiferro-pseudospin (PAF) states show- 
ing the 120° structure of (a r ) in each plane perpendic- 
ular to the rods q = ^fo(l,±l,±l) for S > 1 - 2\q\, 
and (iii) a planar ferro-pseudospin (PF) state charac- 
terized by (& r ) = (cos 0, sin 6, 0) for S < —1/3 and 
5 < 1 — \q\/2 with an arbitrary angle O. Note that 
the degeneracy along the rods q = g(l, ±1, ±l)/\/3 in 
the PAF phase could be lifted by an order-by-disorder 
mechanism which favors the q = order because of 
the higher degeneracy, or by a longer-range interaction 
which is not taken into account in the present paper. 
The nearest-neighbor spin ice (J, q) = (0, 0) and Heisen- 
berg antiferromagnet (5, q) = (1,0), both of which are 
marked in Fig. 7, show no LRO 6 ' 14 down to T = 
but dipolar spin correlations 14 . We stress again that the 
present mean-field treatment becomes inaccurate around 
the nearest-neighbor Heisenberg antiferromagnet. Note 
that the recently studied 120° antiferromagnetic planar 
model 45 corresponds to the limit of S = \q\ —> 00. 

For Pr 3+ ions, the planar components a r represent 
atomic quadrupole moments, as explained in Sec. II B. 
Note that they are defined in local coordinate frames 



through a r = ((<r r • x r ) : (a r • y r )). Then, the phases 
(ii) PAF and (iii) PF are characterized by the fol- 
lowing quadrupole order. Since the local frames sat- 
isfy Eq. (13), the collinear PF state has a noncopla- 
nar ant iferro quadrupole LRO without any translation- 
symmetry breaking. On the other hand, the 2D PAF 
state has a noncollinear alignment of atomic quadrupole 
moments in each [111] layer. It exhibits a noncopla- 
nar ferroquadrupole (FQ) order having a finite uniform 
quadrupole moment pointing to the direction of q for 
q < or a coplanar 120° antiferroquadrupole order for 
q > 0. These can be directly shown by using Eqs. (12) 
and (33) with or without the replacement of n' by n for 
q < or q > 0, respectively. 



IV. SINGLE-TETRAHEDRON ANALYSIS 

It is instructive to investigate the quantum interplay 
of a* and a^ ,y in the model given by Eq. (18) within 
a single tetrahedron. A similar analysis on a model for 
Tb 2 Ti 2 7 16 has been employed 32 . 

In the classical case of f3 = 7 = and thus S = q = 0, 
there appear three energy levels in a single tetrahedron: 

• Sixfold degenerate "2-in, 2-out" configurations | ± 
X), I =LY), and | ± Z) have the energy — 2J n n . and 
are characterized by the direction of the net Ising 
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moment on the tetrahedron T, 
Air = ^2 rrlr = M) ^ 



(35) 



reT 



reT 



which points to ±X, ±V, and ±Z directions in 
the global frame, respectively. Here, Mo is the local 
moment amplitude introduced in Eq. (20). 

• Eightfold degenerate "3-in, 1-out" and "l-in, 3- 
out" configurations have the energy 0. 

• "4-in" and "4-out" configurations |4+) and |4-) 
have the energy 6J n . n . 

With nonzero j3 and 7 and thus nonzero 5 and 
the Hamiltonian Eq. (18) can be diagonalized on a sin- 
gle tetrahedron to yield the following set of eigenvalues 
and eigenstates for a singlet, three doublets, and three 
triplets; 

(i) an Ai g singlet which is a superposition of the six 
"2-in, 2-out" configurations 23,32 ; 

E Alg = -2J n . n .(l-4J), (36a) 
\*A lg ) = ^^(|rX) + |rF) + |rZ)). (36b) 

(ii) two Eg doublets which are superpositions of both 
"2-in, 2-out" and "4-in"/ "4-out" configurations 23 ; 



1*1.) 



E Eg = -2 J„.„. (V (2 + 5) 2 + 6<? 2 - 1 + j) , (37a) 
-| (e^\rX) + e-^\rY) + \rZ)) 

+ c'|4 X ), (37b) 

and 

E' Eg = 2 J n . n . (y(2 + S) 2 + 6q 2 + 1 - (5) , (38a) 



V6 



(e^\ T X) + e-^rY) + \tZ)) - c\4 X ), 



r=± 



(38b) 



with a sign x = ± an d dimensionless functions c' = 
[y/Eq/2] [(2 + d) 2 + 6q 2 + (2 + d) y/(2 + d) 2 + Qq 2 }' 1 
and c = yl — c' 2 , 

(iii) a Ti n triplet described with the antisymmetric su- 
perposition of "2-in, 2-out" configurations 

E Tlu = -2 J n . n ., (39a) 

^=^r(|rX),|rF),|rZ)), (39b) 



and 



(iv) two triplets and a doublet purely comprised of "3- 
in, 1-out" and "l-in, 3-out" configurations, whose 
energy levels are given by — 2J n . n .(<5 ± 2q) and 
6 J n .n. 8, respectively. 



In our case of J n . n . > 0, the ground state is given by 
either a singlet \^A lg ) (Eq. (36b)) or a doublet \^e 9 ) 
(Eq. (37b)), depending on whether 5 is less or greater 
than 



S B (q) = -(Vl+7 2 
respectively, as shown in Fig. 8. 

q 



i)A 



(40) 




2120 singlet 



-1.0 



-0.5 



doublet (2120+41/40) 



0.0 



0.5 



■^6 



1.0 



FIG. 8: (Color online) The symmetry of the ground states of 
the effective Hamiltonian l-L e E in the space of 5 and q in a 
single-tetrahedron analysis. 




FIG. 9: (Color online) (a) Outward normal vectors (green 
arrows) of the surfaces of the tetrahedron, used to define the 
chirality kt- (b) Solid angle subtended by four pseudospins 
cr r . . (c) Distribution of the tetrahedral magnetic moment 
A4t in a cooperative ferroquadrupolar state with (Qt z ) > 0. 
The arrows represent the lattice deformation linearly coupled 
toQF. 

In the above chiral representation of the doubly degen- 
erate Eg states, 1^ ), the index \ — ± represents the 
sign of the net pseudospin chirality of the tetrahedron, 



1 T 

— ^ ^ <T ri • <T r2 X <T r 



through the relation 



(41) 



(42) 
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O ^ O 



-1/3 







1 5 

q = * I o q - R Xi k y , or k? A N.N. spin ice 

R=1, |=+1 or -1, X=+1/-1 * N.N. Heisenberg AF 

l=+l E=+1 



FIG. 10: (Color online) The symmetry properties of the 
ground states obtained with the exact diagonalization of the 
cube in the periodic boundary condition. The dashed curve 
is the boundary between the "2-in, 2-out" singlet and the "2- 
in, 2-out" +"4-in/4-out" doublet for the ground state of the 
model on a single tetrahedron, as shown in Fig. 8. The point 
X = (6, q) = (0.51, 0.89) for Pr 3+ is also shown. 



Here, the summation over the sites t\,t 2 , r 3 on the tetra- 
hedron T is taken so as they appear counterclockwise 
about the outward normal to the plane spanned by the 
three sites for each triangle [Fig. 9 (a)], (kt) gives the 
solid angle subtended by the four pseudospins [Fig. 9 (b)]. 

In this level of approximation, the "2-in, 2-out" and 
"4-in" /" 4-out" configurations are totally decoupled from 
the "3- in, 1-out" /" 1-in, 3-out" , because any single pseu- 
dospin cannot be flipped by the Hamiltonian Eq. (18) 
within a single tetrahedron. This decoupling is an arti- 
fact and a drawback of the single- tetrahedron analysis, 
which is resolved by larger system-size calculations in the 
next section. 



V. NUMERICS ON THE 16-SITE CUBE 

Next, we perform exact-diagonalization calculations of 
the model given by Eq. (18) on the 16-site cube [Fig. 4] 
in the periodic boundary condition. Because of the lack 
of the total pseudospin conservation, this system size al- 
ready gives a large Hilbert space, though we can exploit 
the following symmetry operations; 

(i) the even-odd parity of the total pseudospin, E = 

(ii) the translations T(R) by fee lattice vectors R = 

Si=l,2,3 n *- R * With inte g ers n ii 

(iii) the spatial inversion / about a site, 

(iv) the threefold rotation R about a (111) axis. 



Among these symmetry operations, E, / and R com- 
mute with each other. In general, T(R) commutes with 
I G and E but not with / and R. In our case of the 16- 
site cube with the periodic boundary condition, however, 
there exist only two nonequivalent translations T{Ri) 
with the fee primitive lattice vectors Ri (i = 1,2), since 
f(R 1 )f(R 2 ) = f(R 2 )f(R 1 ) = f(R 3 ) in this case, 
and these two translations eventually commute with /. 
Therefore, we can adopt the following set of commuting 
operators, i7 e ff 7 T(Ri), T(R 2 ), and I. R can also be 
used only in the translationally invariant manifold where 
both T(Ri) and T(Ri) have the eigenvalue 1. 

In Fig. 10, we show the symmetry properties of the 
ground state in the parameter space of 5 and \q\. Special 
points corresponding to the nearest-neighbor spin ice and 
the nearest-neighbor Heisenberg antiferromagnet are de- 
noted as a triangle and a star, respectively, where dipolar 
spin correlations appear without any LRO 6 ' 14 . At finite 
S and/or q, there appear four regions in the parameter 
range 5, \q\ < 1. The boundary S = S B (q) [Eq. (40)] be- 
tween the Ai g singlet and the E g doublet ground states 
in the single-tetrahedron level, which is shown as the 
dashed curve, almost gives one of the boundaries. On 
the left-hand side of the curve, i.e., 8 < #b(<?), the fer- 
roic pseudospin exchange (8 < 0) stabilizes a rotationally 
and translationally invariant singlet ground state with 
the even parity J = E = +1 for both the spatial in- 
version and pseudospin parity (filled blue squares). The 
mean-field result of the collinear ferroic LRO of the pla- 
nar components of pseudospins should be realized when 
8 is negatively large. Therefore, it is plausible to as- 
sign most of this region to the PF (AFQ) phase. The 
ground state having the same symmetry appears in the 
case of antiferroic pseudospin exchange coupling 8 > 
when \q\ is much less than 8. Noting that the U(l) spin 
liquid is stable against a weak antiferroic pseudospin ex- 
change interaction 12 and the mean-field result also gives 
a macroscopically degenerate spin-ice state, this could 
be assigned to a U(l) spin liquid 12 or a quantum spin ice 
without magnetic dipole LRO. It remains open whether 
it is also stable against a weak ferroic pseudospin ex- 
change coupling, namely, whether the U(l) spin liquid 
might appear even in the case of 8 < Sb(q)- 

On the other hand, increasing \q\ out of the above two 
regions changes the ground state from the singlet to sex- 
tets. The sixfold degeneracy of the ground states are 
described by a product of (i) the double degeneracy char- 
acterized by the eigenvalues +1/— 1 for the spatial inver- 
sion / and (ii) the threefold degeneracy characterized by 
the three sets of eigenvalues, (1, —1, —1), ( — 1, 1, —1), and 
(-1,-1,1), for the translation (f(R 1 ),f(R 2 ),f(R 3 )), 
or equivalently, the wavevectors kx = ^(1,0,0), fcy = 
^(0,1,0), and k z = ^(0,0,1). In fact, the singlet- 
sextet transition occurs in two steps. The singlet ground 
state is first replaced by the submanifold of the above six- 
fold degenerate states that has an even pseudospin-parity 
E = +1 (filled red circles). With further increasing |g|, 
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it is replaced by the other submanifold that has an odd 
pseudospin-parity E = — 1 (open red circles). Though 
these states might have an antiferroic LRO of planar com- 
ponents of pseudospins as obtained in the mean-field ap- 
proximation in Sec. Ill, the determination of a possible 
LRO in these regions is nontrivial within calculations on 
a small system size. The particular case of S = 0.51 and 
q = 0.89 which we have found for Pr 3+ is also located 
in the region of the six- fold degenerate ground state, as 
shown in Fig. 10, with the energy ~ — 8.825J n . n . per 
tetrahedron. In the rest of this paper, we will investigate 
magnetic dipole, quadrupole, and chiral correlations in 
this particular case. 



A. Magnetic dipole correlation 

First, we calculate the magnetic dipole correlation, 

r,r i,j ' ' 

(43) 

averaged over the degenerate ground states. For non- 
Kramers ions such as Pr 3+ , this quantity is relevant to 
the neutron-scattering intensity integrated over the low- 
energy region below the crystal-field excitations from the 
atomic ground-state doublet Eq. (7), while for Kramers 
ions, the transverse components &^ ,y must also be taken 
into account. In Fig. 11, we show the profiles of (a) 
S(q)/M§ and (b) S(q)/M$ ■ F Pl s + (\q\) 2 for q = ^(hhl) 
with F Pr 3+(q) being the form factor for Pr 3+ . It exhibits 
maxima at (001) and (003) as well as at (§§0), and a 
minimum at (000), as observed in the dipolar spin ice 9 . 
Note however that the calculated profiles are constructed 
from the on-site, nearest-neighbor, and second- neighbor 
correlations, which gives a good approximation when the 
magnetic dipole correlations remain short-range. Obvi- 
ously, when the spin correlation length is longer, we need 
calculations on a larger system size, in particular, around 
the wavevectors such as (111) and (002), where the pinch- 
point singularity 14 ' 15 appears as in the case of classical 
spin ice 48 . Nevertheless, the failure of the strict "2-in, 
2-out" ice rule can broaden the singularity. 

At the moment, there are not so many experimental re- 
sults on the magnetic dipole correlations in PT2TM2O7. 
The only currently available one is a powder neutron- 
scattering experiment on P^S^Or 19 . It reveals the ab- 
sence of magnetic Bragg peaks and the enhanced low- 
energy short-ranged intensity around \q\ ~ — ~ 0.5 A -1 
and ~ Qz. ~ 1.5 A -1 . These features can be explained in 
terms of the peak positions in our calculated profile which 
are similar to the dipolar spin ice, as shown in Fig. 12. 
This experiment also shows the quasielastic peak width 
~ 0.1 meV saturated at 0.2 K 19 . Such large dynamical 
spin relaxation rate ~ J n . n . can be attributed to the ap- 
preciable quantum nature of the Hamiltonian, i.e., large 
S and q. 




-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 
(hhO) 

FIG. 11: (Color online) Calculated neutron scattering profile 
(a) S(q)/Mg and (b) S(q)F Pr s + (q) 2 /Mg for q = ^(hhl), 
with the form factor F Pr 3+(q). 




0.0 0.5 1.0 1.5 2.0 

q(1/A) 

FIG. 12: (Color online) Poweder neutron-scattering in- 
tensity. Theoretical curves with/ without the form fac- 
tor (blue/magneta curve) and the experimental results on 
Pr 2 Sn 2 0r from Ref. 19. 



B. Magnetization curve 

Next, we show the magnetization curve. The applied 
magnetic field H || [111] partially lifts the ground-state 
degeneracy associated with the inversion symmetry: it 
splits the energies of the 7-odd (-) and J-even(+) ground- 
state manifolds. Then, the ground state has the J-odd(-) 
property. In both submanifolds, the magnetic suscepti- 
bility is finite, as seen from the slope of the magnetization 
curve shown in Fig. 13. This is consistent with a finite 
magnetic susceptibility and thus the negative Tew in 
Pr2Zr207 35 and P^I^Or 21 , and with the absence of an 
internal magnetic field in P^I^Or 37 . The ground-state 
7-odd (-) magnetization curve shows a small step or dip 
around /i#i7/ 'Jn.n. ~ 1-8, in comparison with that of the 
7-even(+) excited state. This indicates that this struc- 
ture develops upon cooling. These results qualitatively 
agree with the experimental results indicating the ab- 
sence and the emergence of the metamagnetic transition 
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FIG. 13: (Color online) The magnetization (left) and energy 
(right) per site calculated for the i-odd(-) ground state and 
the 7-even(+) state by adding the Zeeman term —H • M 
to Heff at the applied field H || [111]. The symmetry of 
the ground-state manifold does not change until a level cross 
occurs at hbH / J n . n . ~ 4.6 to the almost fully polarized state. 
Experimental data on Pr2ir207 at T = 0.5 K and 0.06 K from 
Ref. 21 are also shown for comparison, though they include 
additional contributions from Ir conduction electrons. 



at H c ~ 2.3 T for T = 0.5 K and 0.06 K, respectively, on 
Pr 2 Ir 2 7 21 [Fig. 13]. Requiring that fi B H c /J n ^ - 1.8, 
we can estimate the effective ferromagnetic Ising coupling 
as J e ff ~ J n .n. ~ 0.84 K. However, the magnitude of the 
magnetization is overestimated by about 25 %, proba- 
bly because the experimental results on P^I^Or include 
contributions from Ir conduction electrons. Experiments 
on the single crystals of insulating compounds P^Z^Ot 
are not available at the moment but could directly test 
our theoretical model in the quantitative level. 



C. Multipole correlations 

The sixfold degenerate ground state can be written as 
a linear combination, 



i=X,Y,Z I=± 



(44) 



where c^j being complex constants satisfying the nor- 
malization condition X^/I c ij| 2 = 1- Here, we have 
introduced a ground state associated with ki for 

both I = +1 and J = — 1, which shows a finite cooper- 
ative quadrupole moment defined on each tetrahedron, 
(^i,/|Q J T J >i,/) = 0.0387M 2 ^-, where 




A4, T 8ij, i,j = X, F, Z (global axes). 

(45) 



FIG. 14: (Color online) Upper/lower panels: domi- 
nant/subdominant forms of quadrupole-quadrupole correla- 
tions (QtQt') between the tetrahedrons T and T' dis- 
placed by R = Ri (a), R2 (b), and R3 (c) in the coop- 
erative ferroquadrupolar state with the wavevector kz and 
{Q Z T Z ) / 0. In particular, A Hl ,i = (-0.803,0.274,0.529), 
Ar 2 ,i = (0.274,-0.893,0.529), and Ar 3 ,i = (l,-l,0)/^2. 
Red and blue regions represent the positive and negative val- 
ues of Qr,i. 



Namely, for instance, in the ground state sector associ- 
ated with q = kz, the net magnetic moment A4t m each 
tetrahedron T is inclined to point to the ±Z directions 
with a higher probability than to the ±X and ±F, as 
shown in Fig. 9 (c). This reflects a C3 symmetry break- 
ing in the choice of the ground-state sector q = k^. Thus, 
if we impose the C3 symmetry to the ground state given 
by Eq. (44), it is of course possible to cancel the cooper- 
ative quadrupole moment, (GS\ Qj>\GS) = 0. However, 
it is natural to expect that the discrete C3 symmetry is 
eventually broken in the thermodynamic limit. There- 
fore, we study properties of the ground state having a 
particular wavevector kz and thus a direction for the fi- 
nite cooperative quadrupole moment (GS\Q ZZ \GS) 7^ 0. 
Note that this state shows not only axial alignments of 
magnetic dipoles but also a broken translational symme- 
try, and can then be classified into a magnetic analog of 
a smectic (or crystalline) phase of liquid crystals 39 . 

The cooperative quadrupole moment {Qj) linearly 
couples to a lattice distortion to be verified experimen- 
tally: the four ferromagnetic bonds and the two anti- 
ferromagnetic bonds should be shortened and expanded, 
respectively, leading to a crystal symmetry lowering from 
cubic to tetragonal accompanied by a compression in 
the direction of the ferroquadrupole moment [Fig. 9 (c)]. 
This also indicates that a uniaxial pressure along [100] 
directions can align possible domains of quadrupole mo- 
ments. Experimental clarification of magnetic dipole and 
quadrupole correlations in by NMR is in- 

triguing. 

Next, we perform numerical calculations of equal-time 
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spatial correlations of cooperative quadrupole moments, 
(QtQt>), between the tetrahedrons T and T' displaced 
by R which take matrix form in i and j. To characterize 
the real-space correlations, we diagonalize this matrix to 
obtain the two correlation amplitudes, 

fg„ = (Qt, R ,„Qt',r^, (46) 
having orthogonal forms of quadrupoles, 

Qt, r ,,= J2 x kM ( 47 ) 

i=X,Y,Z 

where = (A^ ?/u , M , A^) with /i = 1, 2 is a set of 
orthonormal vectors satisfying J2i=x y z ^R,fi^R,u = <W 
and ^2 i=x y z ^R /i = 0- Figures 14 (a), (b), and (c) show 
the contour plots of A^ ^X 2 + A^ ^Y 2 + A^ ^Z 2 , which 
represent the diagonalized shapes for the quadrupole- 
quadrupole correlations (Q^Q^,) between the tetrahe- 
drons T and V displaced by R x = (0, a/2, a/2), # 2 = 
(a/2,0, a/2), and = (a/2, a/2,0), respectively. Here, 
red and blue colors represent the positive and negative 
values. The upper and lower panels correspond to the 
forms Qr^i and Qr2 showing the larger and smaller cor- 
relation amplitudes, respectively. There exist dominant 
ferroquadrupolar correlations shown in the upper panels 
of (a) and (b), both of which favors ferroquadrupole mo- 
ments along the z direction. They prevail over antiferro- 
quadrupole correlations shown in (c), and are responsible 
for forming the ferroquadrupole order (QfO 7^ 0. 

To gain insight into a "chiral spin state" observed in 
Pr2lr2 07, we have also performed numerical calculations 
of the chirality-chirality correlation (k,tk>t') between the 
tetrahedrons at T and T' . Note that the chirality kt 
is a pseudospin chirality defined through Eq. (41), and 
is not a simple one defined only with the Ising dipole 
moments a%z r . It turned out that this pseudospin chi- 
rality correlation {kt^t') is weakly ferrochiral between 
the tetrahedrons displaced by Ri and i^, which corre- 
sponds to Figs. 14 (a) and (b). On the other hand, it is 
strongly antiferrochiral between those displaced by R3, 
which corresponds to Fig. 14 (c). Namely, the pseudospin 
chirality, which is a scalar quantity defined on the tetra- 
hedrons forming a diamond lattice, dominantly shows an 
antiferrochiral correlation on the nearest-neighbor pairs 
of the same fee sublattice of the diamond lattice. This 
points to a strong geometrical frustration for a chirality 
ordering. The fate of this pseudospin chirality correlation 
should be examined by further investigations, which may 
open an intriguing possibility of a chiral spin liquid 3 . 

VI. DISCUSSIONS AND SUMMARY 

The effective quantum pseudospin- 1/2 model is quite 
generically applicable to other pyrochlore magnets as- 
sociated with rare-earth magnetic moments, though the 
values of three coupling constants for non-Kramers ions 



and four for Kramers ions may depend largely on the 
materials. In this paper, we have concentrated on novel 
quantum effects in the case of non-Kramers ions, in 
particular, Pr 3+ ions, where we expect the most pro- 
nounced quantum effects among the rare-earth magnetic 
ions available for magnetic pyrochlores 10 ' 33 . The quan- 
tum effects may result in two different scenarios, depend- 
ing on the values of coupling constants: (i) a quantum 
spin ice where the quantum- mechanical mixing of "3-in, 
1-out" and "1-in, 3-out" could be integrated out to bear 
quantum effects in magnetic monopoles, and (ii) a ferro- 
quadrupolar state that replaces the spin ice because of 
a quantum melting. We have obtained ferroquadrupo- 
lar state for the case of Pr 3+ , whose magnetic properties 
explain currently available experimental observations in 
Pr2Sn207 and P^I^Oy. Note that long-distance prop- 
erties are still beyond the scope of our present calcula- 
tions on finite-size systems. Further extensive studies 
from both theoretical and experimental viewpoints are 
required for the full understanding of nontrivial quan- 
tum effects in these systems, in particular, Pr2TM207 
and Tb2TM 2 07. Also from a purely theoretical view- 
point, it will be an intriguing and urgent issue to clarify 
the fate of deconfined magnetic monopoles under the cir- 
cumstance of large quantum- mechanical interactions we 
derived. 

Effects of coupling of localized /-electrons to conduc- 
tion electrons on the transport properties are left for a 
future study. A coupling of Pr moments to the atomic 
and/or delocalized orbital degrees of freedom of conduc- 
tion electrons allows a flip of the pseudospin- 1/2 for the 
magnetic doublet of Pr ion. This could be an origin of 
the resistivity minimum observed in P^I^Or 36 . 
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Appendix A: Crystalline Electric Field 

Here, we give the formal expression for the crystalline 
electric field (CEF) acting on a Pr site. In this Appendix, 
the Pr site is taken as the origin for convenience. The 
coordinate frame is chosen as that defined in Fig. 1 (b). 

We consider the CEF created by nearby ions; two oxy- 
gen ions at the 01 sites, six TM ions, and six oxygen 
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ions at the 02 sites. 



CEF from Ol sites 



Two oxygen ions at the 01 sites are located at ±^^z 
and produce the Coulomb potential 



Uoi(r) = qoi ^2 
16<?oi 



V3a 



E 

£=0 



r=± \r 

' 8r 



>/3- 



2£ 



47T 



-Y&(Sl r ), (Al) 



at a position r from the Pr site, where r = \r\ is assumed 
to be smaller than y/3a/8, Q r represents the spherical 
coordinates of r, and qoi ~ — 2\e\ is an effective charge 
of the oxygen ions at 01 sites. 



CEF from TM sites 



at a position r from the Pr site, where r = \r\ is 
assumed to be smaller than the Pr-02 bond length 

b 2 = ^3(^-^/2 + 3^)^ n 2 = (02,0) with 2 = 

arctan(^(| — 77)), and q<j2 ~ — 2|e| is an effective charge 
of the oxygen ions at 02 sites. 



Matrix elements between the /-electron 
wavefunctions 



To take into account the orbital dependence of the 
CEF for / electrons (7 = 3 and mi = — 3, • • • ,3), it 
is sufficient to include the following terms for U(r) = 
Uoi(r) + U T M(r) + U 2(r); 



Six TM ions are located at |t/ and its symmetry re- 
lated points obtained by successively applying the sixfold 
rotation R 6 about z, and produce the Coulomb potential 



£=1 



Utm(t) = q T M ^ 



1 



n=0 



(r)Yg(n r ) + u^*(r)Y 6 -«(n r ), (A4) 



12q TM ( 2r \ 47r v 6m*ro W^tO \ 

^=0 V 7 6|m|<2^ 
^ 7 m=±l 



where 



(A2) 



at a position r from the Pr site, where r = |r| is assumed 
to be smaller than a/2, Dtm = (f , §), and gxM ~ +4|e| 
is an effective charge of the TM ions. 



3. CEF from 02 sites 

Two oxygen ions at the 01 sites are located at 
=b( V^2( § — v) x + V z ) an d their symmetry-related points 
obtained by successively applying the threefold rotation 
Rs about z, and produce the Coulomb potential 



U02 (r) = QQ2^2 X 



2 

L 



1 



T=± n J \r - r(V2(l - r,)R§x + rjz)\ 



b 2 ^ U02 / 4^ + 1 ^ 

U ^ ^=0 V ^ 3|m|<2^ 



(A3) 



u° 2 £(r) = 



16qoi ( 8r 



2i 



47T 



+ 



V^a VvW V 4£ + l 



12<? TM / 2r 





a 




+ 6g 2 




^02 


4 3 (r) = 


6(702 / 
^02 V 


% ±6 (r) = 


6<7o2 / 
^02 \ 



2i 



47T 



2^ 



2£ 



+ 



Ylq T M 



47T, 

13 



^6 ±6 (^02) 



2r\° 4tt 



F 6 ±6 (^ tm ). (A5c) 



a 7 13 



Then, their nonvanishing matrix elements of 



CEF 



dnY™ i *(n)u(r)Yp(n) (A6) 
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are obtained as 



Tz rn,rn 
CEF 



y0,3 

CEF 



.1 . / 5 ?1 i 3_-0 

r a 2 oo /- "4 



JL /5 S + 1 S 



5 n° 



2 / 5 -0 , 3 „-0 



_ y -3,0 
CEF 



66 


v / 13^ : 


rn 


= ± 




5 


n 




m 


= ± 




25 


22 


v / 13^ : 


m 


= ± 




50 


33 


n/T37t 



(V C 3 E° F )* 

3 /7 



(m = 0) 
(VcV)* 



22 V ;r W ' + nvi^' 



(A7a) 



(A7b) 



y-1,2 _ _ _ v -2,l _ _ ( V h-2\ 

V GEF — ^ V CEF J — V CEF — ^ CEF J 

i FT 



, 5 / 21 _, 
HV2^ 4 -llV^ 6 ' 



2 3 f 3 3 \ * / ^ 6 

^CEF = (^CEF J = 429^6 > 



(A7c) 
(A7d) 



where U (r) and represent the radial averages of U (r) 
and u£(r). 



Appendix B: Representation of the manifold in 
terms of single /-electrons 



\Mj = 4a) 



\Mj = 3a) 
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